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Abstract - The diffraction problem associated with a multilayer material slab re- 
cessed in a perfectly conducting ground plane is formulated and solved via the 
Generalized Scattering Matrix Formulation (GSMF) in conjunction with the dual 
integral equation approach. The multilayer slab is replaced by a surface obeying a 
generalized impedance boundary condition (GIBC) to facilitate the computation 
of the pertinent Wiener Hopf split functions and their zeros. Both E x and H z 
polarizations are considered and a number of scattering patterns are presented, 
some of which are compared to exact results available for a homogeneous recessed 
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1 Introduction 


More sophisticated fabrication techniques and design methodologies have lead 
to a greater utilization of microstrip antenna devices with multilayer substrates. 
When these appear on airborne vehicles (e.g., as part of a flush-mounted antenna 
configuration), their scattering properties become important. 

The scattering of a microstrip configuration may in general be understood 
through a characterization of the scattering substructures which comprise it. In 
this paper we investigate the plane wave diffraction from one such substructure, 
the multilayer semi-infinite slab recessed in a perfectly conducting ground plane 
(figure la). To the authors’ knowledge, this problem has not yet been attempted. 
However, several researchers have considered the special case where the slab is ho- 
mogeneous [ Pathak , 1979; Volakis and Ricoy, 1989; Ricoy and Volakis, 1989a]. A 
common feature of these analyses, though, is a requirement to solve for the fields 
internal to the dielectric, making them cumbersome and impractical to implement 
in the case of the multilayered slab. 

The difficulty in working with the internal fields is alleviated here-in by repre- 
senting the multilayer slab as a surface characterized by a plane wave reflection 
coefficient R ( cos <f> 0 ) (see figure lb), where <f> 0 can be extended through analytic 
continuation in the complex plane. This enables us to carry out the analysis in 
a symbolic manner regardless of the inhomogeneity profile of the layer. Hence, 
although our focus in this paper is the multilayer recessed slab, the derivations 
will be applicable to any vertically inhomogeneous slab. 

The problem herein is formulated via the dual integral equation approach [Clem- 



mow , 1951] in conjunction with the generalized scattering matrix formulation 
(GSMF) [Pace and Mittra, 1964]. The GSMF is applied to the recessed stub struc- 
ture of figure 2a, depicting a perfeclty conducting half plane elevated a distance 6 
above a reflecting surface with a perfectly conducting stub recessed a distance d 
away from the half plane edge. This formulation requires the solution to a number 
of individual subproblems. As illustrated in figures 2b-f, they correspond to the 
problems of direct diffraction, mode coupling, mode reflection, and mode radia- 
tion. Once solutions to each of these subproblems have been obtained via the dual 
integral equation method, they can be combined in accordance with the GSMF 
prescription to yield the diffraction for the original structure in figure la. 

Unfortunately, the conversion of the symbolic solution into one of practical use 
proves to be a formidable task when the reflection coefficient of the grounded slab 
is obtained in its exact form. The fundamental difficulties are related to: 

• the factorization or splitting of the associated Wiener Hopf functions into 
components regular in the upper and lower half complex plane, and 

• the extraction of the complex zeros (i.e., the waveguide modes) associated 
with the split functions. 

The pertinent Wiener Hopf functions cannot be factored analytically and one must 
therefore resort to a numerical scheme (e.g., see Ricoy and Volakis [1989a]). Also, 
in solving for the complex roots of the pertinent split functions, it is necessary to 
employ a search algorithm in the complex plane, a process which is numerically 
intensive. We circumvent these difficulties by replacing the grounded slab by an 
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opaque sheet satisfying a generalized impedance boundary condition (GIBC) [ Karp 
and Kami , 1965]. The GIBCs involve higher order derivatives of the normal surface 
fields and reduce to the standard impedance boundary condition when only the first 
normal derivative of the field is kept. They have been shown to provide accurate 
simulations of thick coatings and in general their accuracy is analogous to the order 
to the condition (equal to the highest derivative retained). Their derivation has 
been discussed in Rojas and Al-hekail [1989], Senior and Volakis [1989], Bernard 
[1987], Volakis and Senior [1989], Volakis and Syed [1989] and Ricoy and Volakis 
[1989b], and amounts to the specification of the constants multiplying the field 
derivatives so that the impedance sheet simulation of the grounded slab yields a 
reflection coefficient that best approximates the corresponding exact one. This 
approximate reflection coefficient is cast as a ratio of polynomials in cos <f> or sin <f>, 
making the determination of the complex poles and zeros of the reflection coefficient 
a simple task. As a result, the required Wiener Hopf factorizations can be obtained 
analytically leading to computationally efficient solutions. 

In the first part of the paper we summarize the GSMF procedure. The dual 
integral equation method is subsequently employed to formulate each of the sub- 
problems and the necessary solutions are obtained for both E z and H z polariza- 
tions. These are given in terms of symbolic split functions which are then evaluated 
for the specific case of a multilayer grounded slab by casting the reflection coeffi- 
cient in a form compatible with a GIBC simulation of the slab [Ricoy and Volakis, 
1989b]. Results are given and the accuracy of the GIBC simulation is examined 
by comparison with known results for homogeneous slabs. 
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2 Description of GSMF Procedure 


In this section, the generalized scattering matrix formulation (GSMF) is ap- 
plied to the geometry given in figure 2a. This consists of a perfectly conducting 
half-plane located a distance S above the grounded slab, with a perfectly conduct- 
ing stub recessed a distance d away from the half-plane edge. To concurrently 
treat both the E z and H z polarizations of incidence, the quantities F z and F x are 
introduced. They are defined as 


F z 

F x 


E„ 

E z polarization, 

Z 0 H z , 

H z polarization. 

Z 0 H X , 

E z polarization, 


H z polarization. 


( 1 ) 

(2) 


and from Maxwell’s equations 


F r = 


jjh dF, 
k dy 


where 


— 1, E z polarization, 


vi = 


( 3 ) 


( 4 ) 


( 1, H z polarization. 

The individual problems to be considered in the GSMF prescription [Pace and 


Mittra , 1964] axe as follows: 


1. Evaluation of the direct diffracted field by the substructure in figure 2b due 
to a plane wave incidence. This field can be expressed as 

F“(<f>,<f> 0 ) = / P dd (cosa,cos(f> 0 -6)e- jkopcot ^da (5) 

Jc 

where P ^ is the spectrum associated with the currents induced on the half 
plane and (/>, <f>) are the usual cylindrical coordinates of the observation point. 



Additionally, C is the complex contour composed of the directed line segments 
[0 — jo o, 0 — jO], [0 — j 0, 7r — jO], [n — j 0, tt + joo] in the complex a plane. 

2. Evaluation of the field coupled into the waveguide due to a plane wave inci- 
dence (figure 2c) as illustrated. We denote the field corresponding to the n th 
coupled mode as 

Fzn (&>) = C n (cos 4> 0 \ 8) e~i knX (6) 

where C n (cos <p 0 ; 6) is usually referred to as the coupling coefficient and k n is 
the propagation constant associated with the n th mode. 

3. Evaluation of the modal field reflected at the stub (figure 2d). This can be 
expressed as r mn e jfcT " r where T mn is the stub reflection coefficient of the n th 
mode to the m th mode. 

4. Evaluation of the reflected field at the waveguide mouth due to the n th mode 
(figure 2e). This can be expressed by R mn (6)e~ jkmX where R mn (<5) is the 
reflection coefficient of the n th mode to the m th mode. 

5. Evaluation of the radiated field attributed to the m th mode incident at the 
waveguide mouth (figure 2f). This field can be expressed as 

(4>) = / Pm (cos a; S) e~^ coe ^da (7) 

J c 

where P m (cos a, <5) is proportional to the spectrum of the currents induced 
on the half plane due to the incident m th mode. 



Accordingly, the scattered field by the recessed stub geometry in figure 2a is 
given by 

F z {<f>, 4> 0 , 8, d) = f [P dd (cos a, cos <£ 0 ; 8) + P mod (cos a, cos <j> 0 \ 8, d)] e _J<:oPCOs(a_0) da ; y > 8 
J c 

( 8 ) 

where P mod (cos a, cos <f > 0 ; 8 , d ) is associated with the presence of the stub and in- 
cludes the contribution of the waveguide modal fields. It can be written in a matrix 
form as [Pace and Mittra, 1964] 

Pmod(cosa,cos<t> 0 ;8,d) = [P m (cos a; 6)] r {[/] - [W mn (d)][r m „][W m n (d)][R mn (<5)]} _1 

■[W mn (d)}[r mn ][W mn (d)][C n ( cos fa 6)] (9) . 

in which the brackets signify column or square matrices depending on whether one 
or two subscripts appear, respectively. In addition, [/] denotes the identity matrix 

and \W mn (d)] is the modal propagation matrix whose elements are given by 

{ e~ jknd , m = n 

( 10 ) 

0, m ^ n. 

To obtain the field scattered by the recessed material slab it is only required to 
set d and <5 to 0 in (8) and (9). In this case, [W mn (d = 0)] becomes the identity 
matrix and [r mn ] reduces to [/] or — [7] for H z and E z polarizations, respectively. 

Thus, P mod becomes 

Pmod (cos a, cos 4> 0 ) d = P mod (cos a, cos 4>o\ 8 = 0, d = 0) 

= [P m (cos O')] {[/] — Ui[72 mn ]} -1 [ViC n (cos <t>oj\ (11) 

where 

P m (cos a) d = P m (cos a; 8 = 0) 
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C n (cos 4>o) d = C n (cos 4> 0 \ 5 = 0) 

Rmn = R mn {8 = 0 ). ( 12 ) 

Hence, the field scattered by the vertically inhomogeneous recessed slab can be 
expressed as 

F: ( <j> , <f>o) = / Prs (cos a, cos e ^ C06 ^da ; y > 0 (13) 

«/ c 

where 

P r3 (cos a, cos <j> 0 ) - P dd (cos a, cos 4> 0 ) + P mod (cos a, cos <f> 0 ) . (14) 

The steepest descent method can then be employed to evaluate (13) and obtain 
the diffracted field. 


3 Plane Wave Diffraction and Mode Coupling 


Consider the plane wave 


pi _ g cos <f>o+y sin 4>o) 


(15) 


F' x = u 1 sin<£ 0 e J '* o(xc ™* o+! ' sin * o > 


(16) 


incident at an angle <j> 0 upon the structure depicted in figure 2b. In the absence of 
the perfectly conducting half-plane, the total fields may be written as (for y >0) 


FT 

= K + f; 

(17) 

ppw 

X 

= k+k 

(18) 


with 


F; = R (cos <j> 0 ) e J M ICO8 & >-y “ n ^ < 0 

F; = -u 1 F(cos^ 0 )sin^ 0 e-’ >o(rco, ^- ! '" n ^> 
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(19) 

( 20 ) 
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where /? (cos 4> 0 ) denotes the plane wave reflection coefficient of the grounded slab 
referred to y — 0. For the general case of a plane wave incident at an angle a on a 
vertically inhomogeneous ground plane coating, R ( A = cos a) may be represented 
as 


R(X) = - 


AjXf - 

A(A) + yr=32£(A)J ’ 


( 21 ) 


where A (A) and B (A) are even functions of A, with any branch cuts in A (A) also 
appearing in B (A) and vice versa. 

The introduction of the perfectly conducting half plane at y = 8 generates an 
additional scattered field component F‘ so that the total fields become 


F z = Ff w + Ff (22) 

F x = Ff w + Ff (23) 


This scattered field is due to induced currents on the perfectly conducting half 

plane, and can thus be represented by an angular spectrum of plane waves. A 

suitable representation is [C7emmou>, 1951] 

f c P (cos a) e^ koSsiaa e~^ koPCOa ^~ a ^da y > 8 

fc Q (cos a) e -i k ° 6t]na + R (cos (x - a)) da 0<y 

(24) 


Ft = 


implying 

— f c Vi sin a P (cos a) koS 8in0, e -, * op 008 ^~°^da 

Ff = <( 

f c v x sin a Q (cos a) e ~ jkoSsma |e -■>*<>*> c °« (*+<») _ / 2 ( cos ( x — a)) c -i*oP«>. (*-<») j da 

(25) 

in which P(cosa) and Q (cos a) are the unknown spectral functions to be deter- 
mined from the boundary conditions at y = 0. These are 


<< 5 , 
y > S, 

0 < y < l 
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(Bl) Continuity of the tangential electric fields at y = 6, — oo < x < oc 

(B2) Continuity of the tangential magnetic fields at y = 8, x < 0 and 

(B3) Vanishing tangential electric field on the perfectly conducting half plane at 


and we note that the boundary condition at the slab surface is implicitly taken 
into account by the representation (24) and (25). 

The application of (B1)-(B3) in conjunction with (24) and (25) is straight- 
forward. It results in the set of equations 


P (A) = -v x Q (A) [l - v x R (A) e--' 2 * o5 ' A ^ y ] . 

/ OO 

Q (A) V 2 (A) e~i k ° xX d\ = 0 ; x < 0, 

■OO 


j°° Q (A) v 3 (A) [l - t>, R (A) e - I ' i - s ' /!r5T ] e-’^dX = 


(26) 

(27) 


x > ((28) 


where we have set A = cos a , A 0 = cos <j> 0 and 

( 1, E z polarization. 


n 2 (A) = 


(A) = 


l/\/l — A 2 , H z polarization. 
l/\/l — A 2 , E z polarization, 


(29) 

(30) 


v 1, H z polarization. 

These are sufficient to obtain a solution for Q (A) and P( A). However, before 
proceeding, it is necessary to rewrite certain terms in the integrands of (27) and 
(28) as products of “upper” and “lower” functions, that is, functions free of poles, 
zeros, and branch cuts in the upper and lower half A planes, respectively. In the 
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W 




process of doing so, we introduce the definitions 


„ 2 (A) V v}(\)v;(\) 

(31) 

v 3 (\) t' v 3 (A) v 3 (A) 

(32) 

1 .. .-liMy/tt def L w (A; 6) U w (A; S) 

1 V ' RWe - L,(X)U,W ’ 

(33) 

in which 


L w (A; 6) U w (A; 8) = A (A) (l + + (A) (l - 

u x e _wyr: (¥4)) 

L, (A) U s (A) = A (A) + Vl - A 2 B (A) . 

(35) 


In these, L s (A), L w (A; <5), (A), v$ (A) are lower functions while U s (A) ,U W (A; 6), 

v t (A), u* (A) denote upper functions. We also note that L W U W is a function 
characteristic to the loaded parallel plate waveguide and its zeros correspond to 
modes in the waveguide. On the other hand, L S U S is a function characteristic to the 
grounded slab with its zeros corresponding to the surface wave modes supported 
by the slab. 

Substituting (31)-(35) into (27) and (28), we have 


/ CO 

Q 0) v 2 (A) v 2 0) e~ lk ° xX d\ = 0 ; x<0, (36) 

-CO 

/I « w •* « v; (A) = 

-v} (A.) 1>3 (A.) fl - a; Q ; x > ((37) 


These coupled dual integral equations can be solved for the unknown spectra by 
examining the analytic properties of the integrands, and the reader is referred to 
Clemmow [1951] and Ricoy and Volakis [1989a] for a more explicit description of 
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this process. From (3(5) and (37) the unknown spectra are determined to be 


3(A) = 


P(\) = 


Lw (A 0 ; S) U, (A) 

2*jU w (\-,6) L S (X 0 ) 

~V\L W (A; 6 ) L w (A 0 ; 6 ) 


v t (A 0 ) u 3 (A.) yfl — A^ 
^■(A)«3 + (A)(A + A 0 ) 

4(A 0 )v 3 (A 0 )^ 
c7(A)u 3 + (A)(A + A 0 ) 


e jkoSy/l^ (3g) 
jko8y/l-\%' ^Q'j 


27rjL,(\)L 3 (A„) 

Recognizing that Pdd( A) = -P(A), we may then set 8 = 0 and substitute for the 
polarization dependent functions to obtain 

L w (A) L w (A 0 ) y/i — A\/l — A 0 


^(A) = 
P«( A) = 


2^jL,{\)L,{\ 0 ) (A + A 0 ) 

Lw (A) Lu , (A 0 ) \/l + A-y/T -t- A 0 


(40) 

(41) 


2njL a (A) L, (A 0 ) (A + A e ) 
where the superscripts e and h refer to the spectra associated with E z and H z 
polarizations, respectively. 

To solve for the field coupled into the region 0 < y < S, x > 0, (21), (33) and 
(38) are substituted into (24) to obtain the integral expression 


Ft = 


f°° ^ (A,) u 3 - (A,) y/l-Xl L w (A 0 ; 6) 

J- oo 2 njvj (A) vt (A) (A + A 0 ) L a (A 0 ) L a (A) U w (A; 6) 

2jA (A) - -(y 1 ~ A I + 2 B (A) cos (k 0 yVT=X*) 


I ft-* / '-''o ^<v 0 y Vi S' J | c * d\ ; 0 < J/(42<5 

This can be evaluated by closing the path of integration via a semi-infinite contour 
in the lower-half A-plane. The sum of the residues of the captured poles then yields 


F: = E c n (A„) | jA(X n ) — + B (A.) cos (k.yj 1 - Aj) 


n=l 


^1-A5 


> e 


— jk 0 x\ n 


( 43 ) 


V 
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where {A n } are the zeros of U w (\-,6), 


C (\) = —2 e J * o5 \/ 1-A ° - \/ 1-A " (^£1 fl ^ (A P ) u 3 (A 0 ) ^1 A^ 

L, (An) L a (A.) ^ (A„; *) (An) ® + (A n ) (A„ + A 0 ) ’ 

(44) 


and 


£C(A n ;£) 


dU w 
d A 


A=A„ 


(45) 


Substituting (29) and (30) into (44) with 8 = 0, we obtain the more explicit forms 
for C„(A 0 ) as 


C' (A.) 
Ci (A.) 


2 Xn, (A 0 ) \/l — A oN /1 — A n 

L, (An) L. (A 0 ) Ul (A„) A 0 + An 

2jjn, (Ap) \/l 4~ Ap-y/1 4- A n 

L> (An) La (Ao) U' w (A„) A 0 + An 


(46) 

(47) 


for E z and H z polarizations, respectively. 


4 Reflection and Launching of a Waveguide Mode 


Consider now then n th waveguide mode field 


F,„ = 


F xn = 


jA(\ n ) 


+ B ( A ") cos { k °yyJ 1 - A«) 


o 3 ^oxXft 


V\A. (A n ) cos (k 0 yyfl - A£) + jviB (A„) yjl - A^ sin {k 0 y\J 1 - A*)] e jk ° xX ' 

(48) 


incident at the waveguide opening (figure 2e). The radiated fields due to this 
excitation may be again represented by (24) - (25). Subsequent application of the 
boundary conditions (B1)-(B3) then yields the dual integral equations (with the 
usual transformation to the A plane) 


/ oo 

<?(A)U f (A) *3 (A) 

-oo 


L.(X)U,W 


0 ; x > 0 


0 < y < 8 


(49) 
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/ oo 7y. ( oi^oXA rl 

.9(J)»?(J) t i(A) e ->‘-V» = _ ; *450) 

2 cos (koSy/l - Ay 


where 


A (A n ) , polarization 


t 5 (A n ) , H z polarization. 

The solution of those proceeds in a manner parallel to the previous case. The 
resulting spectra are determined to be 

Q (A) = 7 >{*n)U,(X)v 3 (An) ( 52 ) 

4tt;'cos (k 0 6yJl^\l)L,(\ n )U w {\]6)v 2 (A)u 2 (A n ) (A) (A + A n ) 

P (A) = v i v 4 L w (\ n )L w (A; <?) (A n ) v 3 (A n ) (53) 

4 vj cos (k 0 S y/l - XI) L a (\ n )L, (A) v, (A) u 2 ~ (A n ) v + (A) (A + A n ) 

with T,, u ± , etc. as defined in (34), (35) and (31), (32). 

Substituting (4), (29), (30), (51) and (33) into (53) and setting (5 = 0 we have 

pe/n _ (A n ) (A) (A n ) \/l — Ay/1 — A n 

' ,/T^aI 47rji, (A) L„ (A n ) A + A n (54) 


Pn( A) = B (A n ) 


MA)M^ \/l +~A\/1 + A n 
4*jL, (A) L„ (A n ) A + A n 


corresponding to the spectra for the E z and H z polarizations, respectively. The 
modal field reflected back into the guide may be computed by substituting (52), 
(21) and (33) into (24) and employing the usual transformation to the A plane to 


obtaii 


■ - j: 


2jA(X) 


-u 4 u 3 (A n ) . Lw (A n ; £) 

4 *j cos (k 0 Syfl - A*)vf (A) u 2 - (A n ) u+ (A) (A + A n ) ’ L, (A) L, (A w ) (A; <5) 

sin (k 0 y y/l - A 2 ) , 1 . 

v _ L i O D / \ A f L. /l V2\ -IKoX A J\ 


T^A2 


+ 22? (A) cos (fc 0 y\/l — A 2 j e * k ° xX d\ ; 0 < y < 8 (56) 


As in the case of coupling, this integral can again be evaluated by closing the path 
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F* 

zm 


* sin (k 0 yJl -A2 ) , \ } 

= E i -I - + B (A m ) cos ( Kyjl - A^ ) i 

m=l \/l - A^ V ' I 


of integration in (56) via a semi-infinite contour in the lower half A-plane to obtain 

-jk 0 x\ m 

(57) 

where R mn are the mode reflection coefficients given by 

R _ V4V3 ( A n) e~ ,fca V 1 ~ A " £,i,(A n ;<$) 

cos {kjjl - A 2)n 2 - (A m ) n 2 - (A„) u 3 + (A m ) (A m + A n ) r, (A m ) (A n ) C; (A m ; 6 ) ' 

(58) 

When <5 is set to zero, this reduces to 

^ (^n) L w (A n ) y/l — A m \/l — A n 


R "" A„, + A„ 


= B(A„) 


-^ur (A n ) \A + A m -y/l + A n 


L a (A m ) L a (A n ) U' w (A m ) A m + A n 
for the E z and H z polarizations, respectively. 


(59) 

(60) 


5 Computation of Spectra for Material Insert in a Per- 
fectly Conducting Ground Plane 


We now have all the necessary components required for constructing the spectra 
P za associated with the multilayer slab recessed in a ground plane as defined in (14). 
Substituting (40), (46), (54), (59) and (11) into (14) we obtain the E z polarization 
result 


pe /\ \ a Lw^ (A) L w (A 0 ) \/l Ay^l A c 

rsl ’ o) ~ 2rjL t (\)L a (\ 0 ) A + A 0 


N N 

i+EE 


v e 

mn 


A + A 0 


^1 1 \ A + A r 


where 


( 61 ) 


f/e ^ (^m) Lyj (A m ) 

" m ~ yi_A^i.(A m )I.(A n )£/i(A n ) 


fl - - KV, 


e 

mn 


(62) 
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(63) 


Kn ={[/]-+ 


For H z polarization, the spectra may be obtained by substituting (41), (47), (55), 
(60) and (11) into (14) to find 


Pr\ (A, A.) = 


— L w (A) L w (A 0 ) \/l + \yf\TT 0 
2irjL s (A) L s (A 0 ) A + A 0 


N N frh 

i + EE — 


A + A 0 


n ~i A a + A n \ A + A r 


where 


vl - B (A m ) £< {K) £ (An) \A + W 1 + A ” v ’. 

V^„= {[/]-[R»n));' n - 


A 

mn 


(64) 

(65) 

( 66 ) 


Expressions (61) and (64) can now be substituted into (13) and the resulting 
integral can be evaluated by the method of steepest descents to yield the far zone 
non-uniform diffracted field 


2tt 


o-ik 0 p 


F s (cos <t>, cos <t>o) ~ \l —e 3 * f4 P ra (cos <f>, cos <f> 0 ) , 

fan yj P 


(67) 


where ( p , <f>) denote the usual cylindrical coordinates. In (67), F 3 and P ra refer 
to E 3 and Pf a in the case of E z incidence and to Z 0 H 3 , P? a for H z polarization. 
Although not apparent, (67) is reciprocal with respect to cos <j> and cos<^ 0 , as it 
should. We also note that P ra (A, A 0 ) is a combination of an inhomogeneous solution 
(direct diffracted term) and a sum of homogeneous solutions (modal contribution). 
It may also be easily shown from the asymptotic behavior of (61) and (64) that 
the homogeneous terms do not affect the edge condition. 
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6 Specialization to the GIBC Representation 


To obtain numerical results, we must first provide expressions for the multi- 
layered grounded slab reflection coefficient (i.e. A (A) and B( A)), as well as the 
associated split functions and corresponding complex roots (waveguide modes). To 
accomplish these tasks in a simple manner, we consider the general GIBC approx- 
imation to R( A). This amounts to setting 

N A 

/1(A) = EMl-A 2 )" 

n=0 

N B 

B( A) = X>n( 1-A 2 ) , (68) 

n=0 


where A n and B n axe constants specific to the multilayered slab and are given in 
[Ricoy and Volakis, 1989b]. Introducing (68) into (21) yields 


R(X) = 


Hn =0 A n (1 - A 2 )" - yT^E£fo£n(l-A 2 ) n 


(69) 


LE&^Ml - A’)" + ^/^^En=!o£»(l - A 2 ) n J 

and note that for a given order of approximation, Na and Nb are finite and in 
general Na = Nb or Na = Nb + 1 , with the order of the condition equal to 


Ns = max ( 2Na , 2Nb + 1) 


(70) 


To evaluate the split functions L w (\-,6) and U W (\\S), (68) is substituted into 

(34) and by setting 6 = 0 we obtain 

2\/l — A 2 Y.n~o B n ( 1 — A 2 )" E z polarization 

2 En=o An (1 — A 2 ) 71 H z polarization. 

The spht functions axe then trivially obtained as 


L w (A) U w (A) = 


(71) 


U w (A) = L w (-A) = < 


yj2B^ B \/l — Afl^=i [\/l -~Tn — A] E z polarization 
flJii [v/l - fn - A] H, polarization, 


( 72 ) 



where Im\J 1 — < 0 and {£„} are the Na zeros of J2n=o ^nX" for polarization 

or Y^,n = o A n x n for H * polarization. It is apparent from (44) and (72) that the 


pertinent waveguide mode propagation constants are given by & 0 \/l — £ n - 
In a similar fashion, we may substitute (68) into (35) to obtain 

Ls(\)U,( A) = J2A n V +£5 n ^/^^ 2n+1 


71 = 0 

N S 


71 = 0 




TV =0 


= *2(. + ^ 

71=1 \ 7n > 


(73) 

(74) 


where 


n is even 


Sn = 


A n f 2 

5(„-i )/ 2 n is odd, 


N S a 


(7n) = {zeros of the polynomial ^ ^7- ( — 1) ; ■y 1 }. 

1=0^0 


(75) 


The factorization of (74) is again trivial upon making use of the well known splitting 
germaine to the impedance half plane problem [ Senior , 1952]. Noting that 

1 


1 + rfy/l- A 2 d ~ 


we have, 


U, (A) = L t (-A) = 


K+ (A -,t})K. (A; t/) 

y/So 


n^i^ + (A;l/ 7 n) 


(76) 


(77) 


where 

A' + (A; ,) = K. (-A; ,) = { 


f A-+ (A; 9 ) 


Re(rj) > 0 


{i 7 ? (a - yi - l/>? 2 ) (A; j?) } ifefr) < !), ^ 
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and 


w 


w 


V 


V 




w 


V? 


K+ (cos a; rj) 


yj [tf „ (3 tt /2 - a - 9) , (tt /2 - a + 0)] 2 
*i(»/2) [l + V5cos(^±*)] [l + V2cos ( M2 z£z£ )| 


i?e (tj > 0) . 

(79) 


In the above Im (yjl — I/77 2 ) < 0, 0 = sin -1 (i7),0 < Re(0), and is the Mal- 
iuzhinets function [Maliuzhinets, 1958], whose evaluation in algebraic form has 
been given in Volakis and Senior [1985]. Whereas the zeros of U w (A) represent the 
waveguide modes, the zeros of L s (A) (which are the poles of KS) correspond to 
the surface waves supported by the material layer. Although not required in this 
analysis, these are easily extracted from (78) and (79). 

The expressions (68) through (79) provide a complete description of a GIBC 
implementation and permit the simulation of any multilayered coating. We remark 
that a unique GIBC modeling of a given coating does not exist; in fact one may 
employ GIBCs of substantially different character to simulate the same configura- 
tion. This point is discussed in the following section and some numerical results 
are provided for illustration purposes. 


7 Numerical Results 


In this section GIBC simulations of various material inserts are presented and 
compared with exact results available for the case of a single layer. Due to its 
greater interest, data is presented only for the H z polarization case. The GIBCs 
employed here-in are given in Ricoy and Volakis [1989b]. These are valid for arbi- 
trary multilayer coatings and are therefore suited for this application. In particular, 
these multilayer GIBCs are synthesized by combining the component- layer GIBCs 


'w 
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in an appropriate manner, pointing to the necessity of understanding single layer 
simulations in order to construct multilayer ones. 

Figure 3 shows the far zone pattern of a single layer insert (e = 2 — y .000 1 , y, = 
1.2, t = .2 A) modeled by various “low contrast” GIBCs (i.e. those GIBCs which 
improve as the layer thickness or index of refraction decreases). To illustrate their 
relative contributions to the far zone pattern, the direct diffraction and modal 
contributions have been isolated in figures 3a and 3b, respectively, with the overall 
result presented in figure 3c. We note that for this low contrast GIBC, an 8 th order 
simulation provides a reasonable approximation to the diffraction pattern. 

In Figure 4, both the thickness and the index of refraction have been increased 
in a low contrast simulation of a (single layer) material insert with e = 3.5 — 
j.OOOl, y = 2.0, and r = .4A. In contrast to the previous figure we now observe that 
a 20 th order simulation is required to obtain a converged result. This degradation 
with increasing index of refraction proves typical of low contrast simulations and 
illustrates the need for other types of GIBCs whose performance improves in this 
range of material parameters. 

In Figure 5 the same material insert corresponding to the data of figure 4 is 
simulated with high contrast GIBC (i.e., a GIBC which improves as the index 
of refraction increases or as the layer thickness decreases). In contrast to the 
data in figure 4, we now observe that only a second order high contrast GIBC 
simulation is required to accurately evaluate the far zone scattering. This difference 
in performance between the low and high contrast GIBCs stems from the type of 
approximation employed in their derivation and the reader is referred to Ricoy 





and Volakis [1989b] for a more in depth discussion. Some insight on the type of 
simulation provided by the low and high contrast GIBCs may be gained through an 
examination of the waveguide modes predicted by the different simulations. These 
are presented in Tables 1 and 2 for the single layer simulations corresponding to 
the data in figures 4 and 5. The exact modes given in the table are generated by 
the equation 



We observe that as the order of the low contrast simulation is increased, the data 
in table 1 reveal that the waveguide modes are “picked up” in a sequential manner 
corresponding to increasing n in (80). On the other hand (see table 2), the high 
contrast GIBCs pick up the n = 2 exact mode immediately and then “branch 
off” to pick up the other modes. The discrepency in pattern convergence between 
figures 4 and 5 clearly suggests that the n = 2 mode is the most significant in 
terms of diffraction (for this particular configuration). We explain this physically 
by noting that the n — 2 mode may be resolved into its constituent rays which 
strike the interface at a characteristic angle (say 9™ od ) which is greater than the 
critical angle 9 C of the material insert. On the other hand, the n = 0 and n = 1 
modes are associated with characteristic angles less than the critical angle. This 
implies that upon coupling into the slab, the n — 2 waveguide mode is partially 
transmitted into free space while the lower order waveguides modes remain bound. 

The above hypothesis may be tested by computing exact solutions in which 
the selection of the included waveguide modes parallels the order in which they 
are picked up depending on whether a low or high contrast GIBC simulation is 
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employed. Figure 6 depicts a high contrast simulation of a single layer having t = 
11— j.OOOl,/* = 7, and r = .4 A. The 10-mode result is a pattern obtained by adding 
in modes sequentially as determined from (80), thus paralleling a low contrast 
mode selection scheme. On the other hand, the single mode result contains the 
contribution of only the n = 7 mode (the mode with 0 m > 0 C ), thus, paralleling the 
high contrast mode selection criteria. This clearly verifies that the most significant 
waveguide modes are those that are “visible”, i.e. those with 9 m > 0 C . 

Finally, figure 7 provides a simulation of a three layer insert composed of two 
high contrast layers (e = 11 —j.0001,p = 7, and r = ,4A and e = 3.5 — j.0001,/1 = 
2.0, and r = ,4A) placed beneath a low contrast layer with e = 2 — j'.OOOl, /i = 1.2 
and r = .2A. These are precisely the layers considered earlier in isolation. One 
might, therefore, expect that the order of the GIBC which provided converged 
results for the single layer simulation will also provide an equally acceptable sim- 
ulation when the slab is part of the multilayer stack. For the case at hand this is 
indeed true, as evidenced by the converged 9 th order result. We also remark that 
the presence of the two high contrast layers enhances the modal contribution to 
the total diffraction when compared with the single layer data given in figures 3b 
and 3c. 

8 Summary 

In summary, the scattering from a vertically inhomogeneous slab recessed in a 
ground plane was obtained through application of the generalized scattering matrix 
technique in conjunction with the dual integral equation approach. The solution 
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was specialized to the case of a multilayered slab simulated with a generalized 
impedance boundary condition (GIBC). Results were given for various single layer 
inserts and it was seen that in the case of materials having sufficiently high index of 
refraction, high contrast GIBC simulations converged more rapidly (with respect 
to the order of the GIBC) and performed better than low contrast simulations. 
Finally, results were presented for a lA-thick lossless three- layer insert containing 
both high and low contrast layers. It was shown that the simulation converged at 
the point predicted by the individual layer simulations, suggesting a method for 
constructing multi-layer simulations. 
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Table Captions 


Table 1. 


Table 2. 


Low contrast approximation to waveguide modes for 
a layer with e = 3.5 — j.0001, n = 2.0, and r = .4A. For each low 
contrast boundary condition, the three numbers of the column headings 
indicate the order of the approximation in thickness r, the order 
of the resulting boundary condition, and the total number of modes 
(see [Ricoy and Volakis, 1989b]). 

High contrast approximation to waveguide modes for 

a layer with e = 3.5 — j.0001,/1 = 2.0, and r = .4A. For each high 

contrast boundary condition, the three numbers of the column headings 

indicate the approximation in the index of refraction k - 1 , the 

order of the resulting boundary condition, and the total number of modes. 


Table 1: Low contrast approximation to waveguide modes for a layer with t = 3.5 — j.OOOl,/* — 
2.0, and r = .4 A. For each low contrast boundary condition, the three numbers of the col- 
umn headings indicate the order of the approximation in thickness r, the order of the resulting 
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Figure Captions 


Figure 1. 
Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


(a) Multilayer slab recessed in PEC groundplane. 

(b) Representation of slab as surface with reflection coefficient R. 

Illustration of recessed stub geometry (a) and associated 
subproblems: (b) direct diffraction, (c) mode coupling, (d) stub reflection, 
(e) mode reflection at waveguide mouth, (f) mode launching. 

H z polarization backscatter echo width for a material 

insert with r = .2 A, 6 = 2— j.0001, n = 1.2 modeled by low 

contrast GIBCs (See Table 1 for an explanation of the legend entries). 

(a) Direct diffraction component, (b) Modal component, (c) Composite. 

H z polarization backscatter echo width for a material 

insert with r = .4 A, e = 3.5 — j.0001, n = 2 modeled by low 

contrast GIBCs (See Table 1 for an explanation of the legend entries). 

(a) Direct diffraction component, (b) Modal component, (c) Composite. 

H z polarization backscatter echo width for a material 
insert with r = .4A, e = 3.5 — j.0001, n = 2 modeled by high 
contrast GIBCs (See Table 2 for an explanation of the legend entries). 

(a) Direct diffraction component, (b) Modal component, (c) Composite. 

H z polarization backscatter echo width for a material 

insert with r = .4 A, e = 11. — j.OOOl, fi = 7 modeled with high 

contrast GIBCs (See Table 2 for an explanation of the legend entries). 

H z polarization backscatter echo width for a three-layer 
material insert with (ri = .4A, t\ = 11. — j'.OOOl, n\ = 7.), 

(r 2 = ,4A, e 2 = 3.5 - j.OOOl, /r 2 = 2), (r 3 = .2A, e 3 = 2 — i-0001, =1.2). 

In the legend entry a, 6, c(d, e), a and b denote the approximation in k~ 1 
of the high contrast layers 1 and 2, c denotes the approximation in r 
of the low contrast layer 3, while d is the order of the composite GIBC 
and e is the total number of modes. 




(a) 



Figure 1. (a) Multilayer slab recessed in PEC groundplane. 

(b) Representation of slab as surface with reflection 
coefficient R. 




suwSbtems- ° , H r ® cessed stub geometry (a) and associated 
. P (k) direct diffraction, (c) mode coupling (d) stub reflection 

(e) mode reflection at waveguide mouth, (f) mode launohfng ° ’ 
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Figure 3. H t -polarization backscatter echo width for a material 
insert with r-.2A, e = 2-j.0001, 1.2 modeled 

by low contrast GIBGs (See Ihble 1 for an explanation of the legend entries), 
(a) Direct diffraction component (b) Modal component (c) Composite. 
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Figure 4. H a - polarization backscatter echo width for a material 
insert with r«.4A, e-3.5-j.0001. m- 2 modeled 

by low contrast GIBCs (See Thble 1 for an explanation of the legend entries), 
(a) Direct diffraction component (b) Modal component (3 Composite. 
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Figure 5. H t - polarization backacatter echo width for a material 
insert with t-.4A, e-3.5- j.0001. u-2 modeled 

by high contrast GIBOi (See liable 2 for an explanation of he legend entries), 
(a) Direct diffraction component (b) Modal component (<y Composite. 
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Figure 6. H # -polarization Backscatter echo width for material 
insert with r=.4X,c=ll-j.0001, modeledwith high contrast 
GIBCs (see Ihble 2 for an explanation of the legend entries). 
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Figure 7. H -polarization Backscatter echo width for a three-layer 
material insert with (r«.4X. t ,-11-1.0001, m, -7 ). 

(t,=. 4X. e =3.5-j.0001. /x,=2). (r *-.2X. e,=fe-J.0(J01. m,= 12)- In 
the legena entry a,b,c(d,e), a and o denote the approximation in 
* ‘ of the high contrast layers 1 and 2, c denotes the 
approximation in r of the low contrast layer 3, while d is the 
order of the composite G1BC and e is the total number of modes. 



